Finite Energy of Black Holes in Massive Gravity 



O 



CD 



> 

o 
o 

rn 
o 



X 



Denis Comelli/ Marco Crisostomi,^''^ Fabrizio Nesti,"* and Luigi Pilo^''^ 

^INFN, Sezione Ferrara, Italy 
^ Universita dell'Aquila, L'Aquila, Italy 
^INFN, Gruppo collegato L'Aquila, Italy 
*ICTP, Trieste, Italy 
(Dated: January 21, 2013) 

In GR the static gravitational potential of a self-gravitating body goes as 1/r at large distances and 
any slower decrease leads to infinity energy. We show that in a class of four-dimensional massive 
gravity theories there exists spherically symmetric solutions with finite total energy, featuring an 
asymptotic behavior slower than 1/r and generically of the form . This suggests that configurations 
with nonstandard asymptotics may well turn out to be physical. The effect is due to an extra field 
coupled only gravitationally, which allows for modifications of the static potential generated by 
matter, while counterbalancing the apparently infinite energy budget. 
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In general relativity, the equivalence principle forbids the 
localization of gravitational energy. Indeed, there exists 
no local geometrical quantity that can be considered as 
a candidate for the energy momentum tensor of grav- 
ity. Thus, all the numerous proposals of local quantities, 
representing the gravitational energy momentum density, 
are noncovariant objects (pseudotensors). Nevertheless, 
it is possible to define a notion of quasilocal gravitational 
energy as an integral on a closed 2-surface of a suitable 
covariant density, see [1] for a review and a full list of 
references. Actually, also in Newtonian gravity the total 
energy E is determined by the static gravitational poten- 
tial $ and is proportional to the flux of V<I> on a distant 
sphere surrounding the system. Clearly, if the static po- 
tential $ decreases at large distances slower than 1/r the 
gravitational energy will be infinite. 

Recently, there has been a renewed interest in trying 
to modify gravity at large distances trough a massive 
deformation of GR (for a recent review see [2]) extend- 
ing at the nonlinear level the seminal work of Fierz and 
Pauli [3]. It was shown that to avoid a number of phe- 
nomenological and theoretical issues [4-8] Lorentz sym- 
metry in the gravitational sector has to be broken [9-12]. 
In this framework, an exact class of black-hole-like solu- 
tions with striking properties were found in [13]. They 
feature a new power-like term, in addition to the stan- 
dard newtonian one, i.e. gtt = 1 ^ 2AI/r + 2S r'^ , where 
S" is a new integration constant. While the energy of 
standard Schwarzschild solution coincides with the mass 
parameter M, when 7 > — 1 one would expect that the 
total energy of the solution is infinite. More in general, 
it is interesting to investigate whether a solution of the 
Einstein equations that falls off at infinity slower than 
1 /r can still have finite energy. 

To address this question, it is important to realize 
that to turn the massless spin two graviton into a mas- 
sive one requires the addition of new degrees of freedom 
(DoF). Following an effective field theory approach, one 
can consider GR plus a minimal set of new fields (usually 
called Stiickelberg fields, borrowing the name from par- 



ticle physics), that provide the needed DoF. The model 
building rules are based on symmetries, without both- 
ering about the dynamical origin of the new DoF. On 
the other hand, if a full theory is provided, the effective 
theory can be viewed as the limit in which all the heavy 
modes are decoupled, leaving only the light ones. This 
point of view is well known in particle physics; for in- 
stance, an effective version of the standard model in the 
broken phase can be safely obtained by sending to in- 
finity the mass of the Higgs particle, retaining only the 
would-be Goldstone bosons and the gauge fields. In the 
case of gravity however, playing naively the same game is 
dangerous because the heavy modes also gravitate, and 
removing them can spoil the energy balance of the full 
system. 

In this letter we address the problem of the total energy 
for a new class of black-hole solutions found in massive 
gravity, showing that the contribution of the heavy modes 
regularizes the total energy, making it finite and positive, 
despite the non- Newtonian behavior of the static gravita- 
tional potential. This fact has profound consequences on 
the interpretations of these solutions. In fact, one may 
accept a slow fall-off of the new term, with — 1 < 7 < 2, 
and thus allow for a non-Newtonian gravitational force, 
still with a finite total energy. 

Gravitational Energy in GR. Let us briefly review 
the notion of quasilocal energy in gravity theories. In the 
past, many proposal of quasilocal energies have been put 
forward, in fact it has now been recognized that there 
is an infinite number of them with reasonable physical 
properties. According to the Noether charge method, in 
a general covariant theory, for each diff one can construct 
a conserved charge that is expressed as an integral of the 
relative generators. For this to construction to work, it 
is necessary that all fields that are varying under a diff 
be dynamical, i.e. that there is no external, fixed back- 
ground. The energy will then be given by the value of the 
Hamiltonian, the generator of the translations in time. 
Following the analysis of Chen and Nester [14], each con- 
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served charge Q(^) is made of a bulk term (integrated 
over a t =const. hypersurface S^) and a surface term 
(integrated over the boundary surface 9St) 



d^xB{0 



(1) 



9St 



For diffeomorphism invariant theories the bulk density 
is vanishing on shell, thus the only contribution to Q(0 
comes from the boundary dY,f 

The surface term has to satisfy few requirements 

• The conserved charge has to be functionally diffcr- 
entiable: the variation 5Q is well defined if all sur- 
face terms vanish. That is the surface terms from 
the variations in the bulk should exactly cancel out 
the variations in the boundary term, when suitable 
boundary conditions are used. 

• The boundary term must be covariant. 

The last requirement is fulfilled by Dirichlet and Neu- 
mann boundary conditions and also requires a reference 
background; often taken to be Minkowski space. Ac- 
tually, as shown by Kijowski [15, 16], the choice of a 
reference background is crucial for a proper definition of 
the symplectic space of GR and it is equivalent to fix- 
ing the boundary conditions or the asymptotics for the 
canonical variables. For a scalar field, both Dirichlet and 
Neumann boundary condition correspond to adiabatic in- 
sulation of the system; indeed, no energy flux across the 
spatial boundary is allowed. Neumann and Dirichlet are 
however the only boundary conditions that allow a co- 
variant definition of the symplectic structure of gravity; 
in addition, in the Dirichlet case, the corresponding en- 
ergy is positive definite [18], and gives at null infinity the 
Bondi mass [17]. 

Using Dirichlet boundary conditions for the metric on 
dYit = St, the quasilocal energy is given by 



E = 



1 



327rG Js, 
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where z"" are local coordinates on St, H'^'^ — \/9 9^"^ a-nd 
r^j, is the Levi-Civita connection in a coordinate basis. 
The difference between a quantity / and its background 
value / is denoted by A/. Thus the quasilocal energy will 
depend on the reference metric g and the timelike vector 
associated with time flow. As matter of fact, the energy 
(2) coincides with the one found by Katz at al. [19] using 
a completely different approach. For our purposes it is 
convenient to consider the following metric 

ds^ = -J{r)dt^+tC{r)dr'^+2V{r)dtdr + Br'^ d{?. (3) 

We take as reference the flat background 



with J' and B constants. Thus, if the time flow vector is 



£. — i^i''^) St^ energy is given by 



E = 



Br 
2G 



J 



J 



(5) 



where wc have set F = JfC+T)^ and J- is the correspond- 
ing quantity for the background metric. For the standard 
Schwarzschild metric J = 1//C = 1 - 2 GM /r,V = and 
J = B = 1. When we set = 1 (the time flow coin- 
cides with the Killing vector) the energy coincides with 
the mass of the black hole M . 

Massive Gravity and Exact BH Solutions. An el- 
egant way to provide the needed new DoF for the mas- 
sive phase of gravity, without breaking diff invariance, is 
to work with an extra tensor g^i,, with a kinetic term 
similar to GR, this bring us into the realm of bigravity 
theories whose study started in the '60 (see [20] for early 
references). Once we couple nonderivatively the extra 
spin two to the standard metric field we can explore both 
Lorentz invariant and Lorentz breaking phases of massive 
gravity, working with consistent and dynamically deter- 
mined backgrounds. The action for the modified gravity 
model we will consider is 



S= d^x ^Mi n 



~gMlin-A{gcif/^V , (6) 



where the interaction potential is a scalar function of 
the tensor ~ g^'^gav Clearly the action is invariant 
under diffs. 

For a large family of potentials the equations of motion 
admit the following spherically symmetric solutions [13] 



(where M'^ = IQttG and k = M^i/M^i): 



• For the metric g^j^ 

ds^ = -J{r) dt^ + J"^(r) dr^ + dil^ 



J{r) = l-(-2G 



M 



+ Sr'< 



(7) 



For the metric 



ds 



-J{r) dt^ + k{r) dr^ + 2D{r) dtdr + B r'^dVt'^ 



J{r) = c^uj^ 



B =uj\ 



G 



cM 
r 



K{r) = {c^u:^ + D''{r))/J{r). 



(8) 



where c, uj and 7 are coefficients depending on the 
potential V . In these solutions, besides the standard 
term M/r, there is a power-like modification with a 
new integration constant S which sets its size.^ When 
7 < 2, both metrics describe asymptotically fiat spaces 
iji-^ivpa ^ r^'"'^), namely Minkowski space {rj^^) for g^i, 



^ When 7 = —1, the modification is logarithmic, — ^ log r/r. 
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and the metric fj^,, = a;2Diag(— c^, 1, 1, 1) for g^y. No- 
tice that when c 7^ 1, asymptotically, Lorentz symmetry 
is broken down to rotations [11]. The solution for g^^^ is 
quite similar to the one for (7^1,, however an off-diagonal 
t — r contribution D is present. Its explicit form is not 
essential, but it is important to recall that at large dis- 
tances its behavior is D ~ l/v^ (for 7 < —1). This 
is slower than the typical newtonian weak-field behavior 
and is analogous to the Schwarzschild solution written 
in Painlcve coordinates. Contrary to that case, here this 
term cannot be gauged away; in fact, by the change of 
coordinates dT = dt — D / J dr, one can bring 5^ 1^ to di- 
agonal form, but then g^^ is non diagonal. Incidentally, 
the asymptotic behavior of D implies that the ADM en- 
ergy found in the standard Hamiltonian formulation of 
GR [21, 22] is not well defined [23]. 

A similar solution for g^^ exists also in the Stiickelberg 
effective approach [25, 26]; this time g^i, is replaced by 
a set of four scalar fields. In the decoupling limit of g^j^ 
(k — >■ 00) the fluctuations of the auxiliary metric are 
decoupled, and g^i, is practically frozen to its background 
value fy^jy, modulo gauge. From (7), (8) one can see that 
in this limit, the new term proportional to S survives in 
g^^, while g^^ becomes flat [13]. A preliminary study of 
the thermodynamics of these black hole solutions can be 
found in [24]. 

TotaJ Energy. The computation of the total energy of 
the system described by (6) follows the same steps as in 
GR. The crucial point is that the interaction part in the 
action does not provide new boundary terms in the total 
Hamiltonian (due to the fact that the two sectors interact 
without derivative terms). As a result, the diffs Noether 
charge can be written, on-shell, as the sum 



t(0 



d^3 



St 



(9) 



where B{^) and i5(^) are the single contributions from 
the kinetic action of g and g respectively. The energy of 
the solutions (7) and (8) can be worked out by using the 
general formula (5) applied to g and g. For g we take 
as reference space the Minkowski space 77^1/, while for g 
the natural reference space is its diagonal vacuum metric 
77. Taking as usual ^ = the contributions to the total 
energy coming from g and g are 



E^—il-J)=M-Sr''+' 



E 



2cG 



,2, ,2 



(10) 

.Pj = M + S r''+^ . (11) 



The term 5* 7-'''+^ in both E and E would lead to a diver- 
gent energy in the limit r — > 00 when — 1 < 7 < 2, as 
in the newtonian case. Amazingly instead, the true total 
energy is finite: the power-like terms proportional to S 
cancel exactly, leading to 



Etot = E + E = M + Mc^ 



(12) 



The system has thus finite energy, given just by M and 
M from the Newtonian terms. 



Let us observe at this point that in the effective-theory 
approach the picture is drastically different. As we men- 
tioned, taking first the limit k — > 00 in the solution (7) 
(8), the auxiliary metric gets frozen to its fiat reference 
background 77, while the physical metric g is unmodi- 
fied. Therefore one would attribute zero energy to the 
auxiliary metric and now the total gravitational energy 
appears to be divergent when 7 > — 1. In the decou- 
pling limit thereore, an infinite amount of energy is ne- 
glected, spoiling completely the nice cancellation mecha- 
nism. This also shows that the computation of the energy 
and the decoupling limit do not commute. Clearly, the 
cancellation can take place only if the auxiliary metric is 
dynamic. 2 

It is also interesting to study what happens when the 
total energy is expressed in terms of the body internal 
structure. For the sake of simplicity we considered a 
body of radius Rq, made of constant density fluid [13], 
with total mass m = ^ttRq p. The integration constants 
M, M and S can be determined by matching equations 
(7), (8), representing the external part of the solution 
(valid for r > Rq where Rq is the body's radius) with 
the interior one. The result of the matching is 



and 



M = m + AM, 



M = - ■ 



AM 



AM 



t^grav i 



(13) 



(14) 



with b a suitable constant not relevant for the present 
discussion (see [13]). The mass parameter fi^^^^^ is deter- 
mined by the potential V: generically, it has indefinite 
sign and is of the same order of the graviton mass. The 
nonzero value of S means that the modification of the 
potential is turned on by the matter source m, and the 
appearance of AM in (13) shows that the body mass m 
in the Newtonian term is renormalized, leading to screen- 
ing/antiscreening of mass. This happens because of the 
gravitational interaction between g and g. 

Nevertheless, evaluating eq. (12), we get that the total 
energy coincides with the "bare" mass tti. 



Etot ^ M + Mc^ = m. 



(15) 



This result is striking, since the total energy is indepen- 
dent not only from S but also form the mass renormal- 
ization AM, suggesting that some sort of Gauss Law is 
still valid for the total energy. 



^ We also checked that, for the less physical Neumann boundary 
conditions [11], where SF"^^ = at the boundary, the corre- 
sponding energy -E^f" is in general not finite anymore. When 
7 < —1 one still recovers the finite total energy (12) at large r, 



however Ej^ 
ular E^-- 
< 7 < 2. 



turns out to be divergent if 7 > 



rT+l for -1 < 7 < and E^J" 



-1. In partic- 
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Conclusions. We studied the total gravitational energy 
of the exact spherically symmetric solutions, in a class 
of massive gravity theories where the metric field couples 
in a nonderivative way to an additional tensor field. In 
such solutions the gravitational potential features an in- 
teresting modification $ ~ M/r + Sr^ , which describes 
an asymptotically flat spacctime for 7 < 2. Clearly, when 
— 1 < 7 < 2, the asymptotics is peculiar and in conflict 
with the standard behavior $ ~ l/r for r — )• 00, ex- 
pected in a system with finite energy. Nevertheless, we 
have shown that the new field responsible for the modifi- 
cation also contributes to the total energy that turns out 
to be finite and positive. The point is that the seemingly 
infinite amount of energy due to the r"' term is cancelled 
by a similar contribution from the backreaction of the 
additional field. 

The new field is coupled to gravity only, and while it 
will affect the propagation of gravitational waves [11], it 
cannot be probed by standard matter experiments. The 
standard observers which couple to metric g^j^ will not 
experience any violation of the weak equivalence prin- 
ciple, but they will observe the peculiar behavior of the 
gravitational potential $ at large distances. Just by look- 
ing at the asymptotics of they would be baffled by its 
behavior and led to conclude that the "total" gravita- 
tional energy of the system is infinite, while in reality it 
is not. The total energy is not only finite but also does 
not depend on the constant S appearing in the solution. 
Though the BH-likc solution has hairs (the values of S) 
they have no effect on the total energy. 



When the solution describes the exterior part of a self- 
gravitating body, e.g. a realistic star, the integration con- 
stants M, Af and S can be computed in terms of the 
body's density and radius. While both Af and M show 
in general (anti)screening of mass, and together with S 
arc in general complicated functions of the body shape 
(size), the cancellation mechanism is again operative, and 
as a result, the total energy content of a stellar object is 
just its "bare" mass. 

The picture that emerges is that in modified grav- 
ity theories the gravitational potential can greatly differ 
from our naive expectations. The potential can have a 
softer than expected fall-off, without clashing with the 
basic principle that any finite physical system must have 
finite energy. This we think is remarkable and can change 
our view on gravity especially if the cancellation mecha- 
nism is a general feature of modified gravity. 

It will also be very interesting to explore the cosmolog- 
ical implications of the scenario described in this work, 
and we leave it for a future investigation. Clearly, an im- 
portant open question is the stability of our solutions, in 
particular when the term r'^ dominates at large distances 
(—1 < 7 < 2). The answer is not an easy one, entailing 
the nontrivial study of perturbations around the solution 
(7)- (8), and clearly deserves further study. 
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